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The number of bound states resulting from inhomogeneities in a BCS superconductor is usually 
established either by variational means or via exact solutions of particularly simple, symmetric 
perturbations. Here, we propose estimating sub-gap states using the Birman-Schwinger principle. 
We derive upper bounds on the number of sub-gap states for small normal regions and find a 
suitable Cwikel-Lieb-Rozenblum inequality. We also estimate the number of such states for large 
normal regions using high dimensional generalizations of the Szego theorem. The method works 
equally well for local inhomogeneities of the order parameter and for external potentials. 



I. INTRODUCTION 

It is well known that inhomogeneities in a supercon- 
ducting system may give rise to Andreev bound states. 
Such bound states form in superconductor-normal metal- 
superconductor (SNS) junctions and can affect the trans- 
mission and noise properties of the junction. The pres- 
ence of such states in the cores of vortices is of much 
interest as well. Vortices in s-wave superconductors 
may have Caroli-de-Gennes-Matricon (CdGM) states pQ, 
bound states in vortices in p + ip superfluids where stud- 
ied e. g. in [5]. Zero modes in the core of half-quantum 
vortices in p-wave superconductors [21 H] are of extensive 
recent interest. Such modes are described in terms of un- 
paired Majorana fermions and have non-abelian braiding 
statistics E]. Other inhomogeneities in the system 
have been of interest as well, for example, stripes and 
other structures in the pairing function of a supercon- 
ductor (See e.g. 015]). 

Often, the bound sates are either directly obtained for 
systems of high symmetry, such as a circular vortex, or 
idealized planar SNS junctions or studied using pertur- 
bation theory. In addition, multiple scattering and semi- 
classical (WKB) methods based on the Andreev approx- 
imation have been developed [5J. 

In this paper, we introduce an alternative view on An- 
dreev states, complementing the above-mentioned meth- 
ods. Our approach does not assume a particular sym- 
metry, and works even when the Andreev approximation 
is inadequate (see e.g. 10J) . We base our approach on 
a simple adaptation of the Birman-Schwinger counting 
argument to the superconducting scenario. The Birman- 
Schwinger method has been widely used to study bound 
states of the Shcrodinger equation with a potential. The 
Birman-Schwinger bound has been introduced by Birman 
[TT] and Schwinger [12] as a refinement of previous esti- 
mates by Bargmann |13j on the number of bound states 
in a potential. For a discussion and review see [13]. Re- 
cently the Birman-Schwinger approach was used in the 
context of estimating the critical temperature of BCS 
models in [I5HI7]. 

In this paper we establish the following results: 

1. A bound on the number of sub-gap states in a su- 
perconductor containing a normal region. The bound on 
the number of states below a given sub-gap energy E is 



expressed in eq. (Ill, and holds for general values of the 
chemical potential (i, the gap Ao- For example, in the 
limit of large chemical potential, the inequality is given 
in the simple form 



N E < 



y/jiAfooliA) A 2 , (Ag - 2E 2 ) vol(A) 



^(A§ -£2)3/2 87rM 3/ 2v / A 2 _ E 2 



(1) 



where vol(A) is the volume of the normal region. Clearly, 
the inequality is most effective for energies not very close 
to the top of the gap. 

2. In studying the bound states of Schordinger oper- 
ators with deep potential wells, often Birman-Schwinger 
inequalities yield poor upper bounds to the number of 
bound states. In such cases the Cwikel-Lieb-Rozenblum 
inequality [T8H2"U] may be used to obtain upper bounds 
with power low behavior in-line with the semi-classical 
expectations for deep wells. We establish the Cwikel- 



Lieb-Rozenblum inequality (20) for the Bogolubov de 
Gennes equation. This bound requires more work and 
has the advantage of retaining a simple form for a spa- 
tially varying A, however, for normal regions, and in the 
typical situation of gap energy being small compared to 
Fermi energy, it seems to give a worse upper-bound to 
the number of states than the simple Birman-Schwinger 
inequality. 

3. Finally, using semiclassical tools and the relation to 
Szego theory we supply an asymptotic expression for the 
number of states in large regions and a conjecture (28) 
for the scaling l d ~ l log(kpl) of the sub-leading term. 

The paper is organized as follows: In section [TTJ we 
derive a Birman-Schwinger operator which is suitable for 
the treatment of mean field BCS superconductors. In sec- 
tion lllll we derive the basic bound on the number of states 
in a normal region of a superconductor. We proceed to 
present the Cwikel-Lieb-Rozenblum inequality in section 
|IV[ the derivation of which is technical and deferred to 
the appendices. In section [V] we consider the bound- 
ary correction to the number of states for asymptotically 
large normal regions and explain their connection to the 
theory of Toeplitz operators and Szego theorems. 



II. THE BIRMAN-SCHWINGER PRINCIPLE 
FOR BCS 

To employ this method, we consider a superconductor 
described by the mean field BCS Hamiltonian: 

Hbcs = £„/ 1>t(x)(u(x) - M - ^V 2 )^(z)da; +(2) 
JiA^ix^'^Kx^lix') + h.c.]dxdx' , 

where A is the superconducting order parameter, (i is 
the chemical potential u(x) is a local potential and m 
the effective electron mass. A is related to the anoma- 
lous Green's function through a self consistency condi- 
tion. Here we will assume that u(x) and A are given, 
and represent the effective values of these parameters, 
which are, in principle, measurable directly. 

To diagonalize the BCS Hamiltonian we use the Bo- 
golubov de Gennes (BdG) equation 



H 



BdG 



4>e 
4>h 



= E 






. tt _ , h A 

, tlBdG — I ^* _fr 



(3) 



where ip e , iph are the electron and hole parts of the quasi- 
particle wave function, and E is the quasi particle excita- 
tion energy measured relative to the Fermi energy. The 
solutions of this equation come in +E, —E pairs (un- 
paired states may exist at E = 0). We will mostly con- 
sider h = — V 2 — u + u(x) (Throughout, we take £- = 1). 

For a translationally invariant A and u(x) = 0, the 
spectrum of the BdG operator is continuous, and con- 
sists of the positive and negative energy bands. In an 
s-wave superconductor these bands are separated from 
zero by the gap |A|. The method will also work for 
other gap functions, for example in a p-wave supercon- 
ductor or superfluid with gap A(p) ~ Aq(j> x + ip y ) (as in 
the Anderson-Brinkman-Morel phase of 3 He). Andreev 
states in d-wave superconductors have also been of much 
interest, however in in d-wave superconductors the gap 
function vanishes on the nodal lines and one needs to 
differentiate between resonant states and the rest of the 
spectrum. 

To study the introduction of bound states into the sys- 
tem from an inhomogeneous A (and possibly also u), we 
write the perturbed HsdG as H + W where H is trans- 
lationally invariant. Here W represents the perturbation 
(inhomogenous A, disorder potential etc.). We choose 
a decomposition of the perturbation W into a product 
W±W 2 , and define the Birman-Schwinger operator: 



K E = W 2 



1 



E-H 



W x 



(4) 



The key property of Ke is that it has an eigenvalue \ 
for each eigenvalue E of (Hq + XW). Indeed: 



(H + \W)<p = Ecf> =4> - ( /> = -(Ho-E)- l WiW 2 < l 
A 



(5) 



(W 2 (f>) = W 2 (H Q - E)- 1 W 1 (W 2 <I>) = K E (W 2 ^ 



Thus, W 2 (j) is an eigenvector of Ke with eigenvalue j-. 
The two major differences between the BCS system and 
the Birman-Schwinger kernel Ke typically used to study 
bound states in Schrodinger operators is that Hb<ig is 
not bounded below, and that Ke is not Hermitian, but 
the ideas have been used to estimate the number of states 
within a spectral gap (see e.g. [21]). 

We now invoke the Birman-Schwinger argument: 
Let N EljE2 (V) = number of eigenvalues of H Q + V in 
the interval (Ei,E 2 ), where Ei,E 2 are in the gap, i.e. 
I .Ei |,|-E2 1 < mm p(|A(p)|). Assuming that the energy of 
eigenvalues of Hq + XV depends continuously on A as we 
increase A from to 1, we see that each eigenvalue for 
A = 1 must have crossed the E\ or E 2 line for some value 
of A e (0, 1), as illustrated in Fig. [I] Thus: 

N(E) < number of crossings fot X G (0, 1) = 

^eigenvalues of Ke larger than 1 (6) 
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FIG. 1: The Birman-Schwinger principle: continuity of the 
bound state energy as function of A implies that each bound 
state at A = 1 with energy in (E, —E) corresponds to at least 
one crossing of the E or —E line for some A £ (0, 1). 

Note that the inequality sign in (Tsl) is due to the fact 
that in all cases we consider perturbations that are not 
strictly positive or negative. The curves in Fig. [TJmay be 
non monotonous functions of A and an in-gap eigenvalue 
at A = 1 may correspond to several crossings. Moreover, 
in some situations no bound states will be created even 
for very large A. 

It is important to note that for unpaired zero modes of 
an unperturbed system such as the p x +ip y superconduc- 
tor (for example in the presence of vortex or boundary) , 
the principle will still work. Such a zero mode is stable 
as long as it is far from other vortices, and will not shift 
when adding impurity potentials. 

We now find the explicit form of the Birman-Schwinger 
kernel useful for the study of Andreev bound states in 
a normal domain of general shape embedded in the su- 
perconductor. First, we decompose: A = A + <5A, 
S A = SAiSA 2 , and note that 



SA 
6A* 



£Ai 
SA* 2 



5A 2 
5At 



(7) 



To further simplify the system we assume an s wave 
superconductor. This will be described by A = Ao out- 
side the region and A = inside (Fig. [5]). Thus we write 
<5Ai = (j x xa{x) and 8A2 = — AoIxyt(#)j where I is the 
identity matrix, and xa(x) is the characteristic function 
of the region A, so xa(x) = 1 if x € A and other- 
wise. The Birman-Schwinger operator Q can be written 
explicitly. Assuming real Aq for simplicity, we find: 



An 



VAo 2 -^ 2 



Ka,e(x,x') = 
(a x E - Ao)xA(z)Img + (x - x')xa(x') 



-iA Q a v XA(x)Re g+ (x - x')xa (x') 



where 



g±(x-x') 



~±ik±\x~x'\ 

Ait\x — x'\ 



(8) 



(9) 



where k± = sj /j, ± i\/A 2 — E 2 , and Re k > 0. 

Note that from the numerical point of view studying 
the Birman-Schwinger kernel & has an immediate ad- 
vantage: the kernel is only supported on the region A. 
Thus, computation of the spectrum of bound states re- 
quires merely finding the eigenvalues of this operator re- 
stricted to the region A. Since the operator in ^ is 
known analytically, one can numerically compute eigen- 
values by simply discretizing the region A. 

Let us proceed and consider limits where the applica- 
tion of these ideas is particularly interesting: 



III. NORMAL REGIONS: 

Here we use Ka,e to obtain an upper bound to the 
number of states. The basic inequality is: 

Ne(V) < Ys\\ l i n (\)=E 1 or E 2 for some AG (0,1) * — ( 1( ^) 



< 



Z^A|^„(A)=_Ei or E 2 for some AG (0,1) A 2 

V — + V 

t-~d A eigenvalue of Kei A 2 £-~< A eigenvalue of Ke2 A 2 

TrK* E1 K E1 + TtK* E2 K E2 



A< 




A=Aq 

FIG. 2: A normal region A inside the superconductor. 



The last inequality is due to the fact that Ke, is, as 
opposed to the usual Birman-Schwinger case, not self- 
adjoint, therefore we used that for a compact operator 

E el genvalue S \ X \ P < Y, singular values \ S \ F for ? > ] (SdlUT- 

Lalesco-Weyl theorem). 

This allows us to write a simple inequality for the num- 
ber of sub-gap states. Plugging Ke,a from ^ in ( 10 ) we 
find: 

N(E)<4AlJ AxA dxdx>(0^\Irng_( X -x')\2 + 

\Reg + (x-x')\A (11) 



As E — > the number of states Ny(E) goes to zero, 
thus, the estimate (11 1 is effective only when the expres- 



sion on the right hand side of ( 11 1 is less than 1 in that 



limit. Let us estimate the size of the region A where this 
holds, assuming a typical gap much smaller than kp. Us- 



ing that for E — > , k± ~ ^/JI ± i^= , we see that 



N(E) < ^ J AxA dxdx 



\x— x'\ 2 



(a^^n/mN-^'D + i) 



(12) 



and so the decay is set by coherence length £ = -^ . This 
means that we need the condition ^ ^/JLAqV ol{A) < 1 to 
make (11) effective at low energies. 



Much like for the Schrodinger equation in 3d, we see 
that for a small enough normal region, the right hand side 
of (111 will be less than one, and so no sub-gap states 



will be present below this energy E. 

Numerically, we have checked the utility of this ap- 
proximation for simple situations on a lattice compared 
to the exact result. The agreement is fairly good for E 
small E < 0.25A . 

Let us estimate the leading volume dependence of the 
bound (fill. Writing 



9±(r) 



3 ±ik±\r\ 

47r|r| 



- ar±ib\r\ 

47r|r| 



(13) 



with: 



: ( M 2 + A 2 - £ 2 ) 1/4 sin[i arctan ^EE] 
(M 2 + A 2 - £ 2 ) 1/4 cos[| arctan V^RE] 



(14) 



By ([11] we have: N(E) < Ai where 

' A 2 n +E~ 



N 1 =4AU AxA dxdx>(^y\Img^-x')\ 2 + 

|Reg + (x-x')| 2N ) (15) 



now, we change to variables x,x' — > x,r = (x' — x). 
Taking the r integral over all space we get that Ai < A{ 



with: 



N[ 



167rA^(^)/ °°r 2 dr(|^|Img_(r)| 2 + |Reg+(r)| : 



- A 2 



vol(A)f-dr^(0^. 



a 2_ E i siir(6r) + cos 2 (&r) 
This last integral can be carried out analytically yielding 



N[ = §- a vol(A){\ + 



-y) 



Finally, we can write N[ as: 

vol(A)A 2 ((2A 2 -E 2 ) v /p, 2 +A 2 -E 2 +p,E 2 



N{ 



47v(A 2 -E 2 )(^ 2 +A 2 -E 2 )- 



-iV^F 



(16) 



(17) 



The upper estimate N[ in ( 16 ) is in fact a good apprx- 



imation to Ni whenever the size of A is larger than 1/a 
the decay scale for the Green's function. 

The most physically common situation for BCS is 
jj, >> Aq. In this limit we have: 



N[ 



JJiA*vol(A) Al(Al-2E 2 )vol{A) 

7t(A 2 -E 2 ) S / 2 + 87T M 3/2^ A 2_ B 2 



(18) 



Taking \x fixed and Ao large, while keeping a — Aq/E 
fixed, we get 



N[ 



2s/2-n{l-u 2 fl- 



(a 2 +2)fi-/A^ 
4\/27r(l-a 2 ) 7 / 4 



(19) 



The bounds (201 and (21) are closely related to the 



bound of [22] on the number of bound states of perturbed 
Dirac operators. The main technical difference with [S3] 
is the presence of a chemical potential and thus a presence 
of a non vanishing Fermi surface [351 . 

Note that the denominator (|A | 2 — E 2 ) 2 in (20) shows 
that this bound is only useful for E < |A |. However, 
most of our interest here is actually with states that are 
not too close to the top of the gap itself. The reason 
for this is that in practice, if the BCS self consistency 
is taken into account, the first states that are affected 
are those with energy very close to the gap. Such states 
will often be washed out from the spectrum, since the 
gap function may decrease slightly in the vicinity of the 
normal region. 

In the case of a system defined by a normal region we 
Ai(:e) = -AqXa(x), we have N(E) < N C lr with : 



N, 



CLR 



= C 



2m 1/2 |A | 5 
|A | 2 -£ 2 ) : 



;VOl(A) 



(22) 



Let us compare the simple estimates CLR and BS esti- 
mates. Consider the large jjl limit (22) vs (18). We have: 



N'JNc 



Ao 



<1 



(23) 



We conclude that for the typical BCS situation, the 
simple Birman-Schiwnger bound ( 18 ) is more effective in 
dealing with normal regions. 



IV. CWIKEL-LIEB-ROZENBLUM INEUQALITY 



Often, the estimate ( 10 ) may be quite poor because 



of large contributions from small eigenvalues of Ke.a- 
This behavior may sometimes be remedied by considering 
other inequalities, such as the celebrated Cwikel-Lieb- 
Rozenblum inequality [T%H2"0"] . 

To do so we consider the BdG equation with a gap 
function A + Ai(x). Using the ideas in [T51 |2"5] and 
adapting to the superconducting scenario we find that 
the inequality 



N(0,E) <C- 



lt 



l/2| An |2 



I, 



\A Q \ 2 -E 2 ) 2 
dx 



|AiW|>i 



is valid for /i > 



Ao /2 |Ai(*)l 3/2 



|Ai(x)| 3 . (20) 



18\/2 



with C a constant of order one. 
The details of this calculation are given in the appendices. 



For the case [i < 



ygFg 



to 



18\/2 



, the inequality is changed 



N(0,E) < 



2 -3/4 C 



/ 



3(A 2 - £ 2 ) 3 / 4 



d.r 



IAxWI^J 



> 3 / 2 l 



AHAiWr + IA^r . (21) 



V. ASYMPTOTIC REGIME, SZEGO THEOREM 

Here we are interested in the large I behavior of the 
eigenvalues of Ke,ia, where I A is the region A rescaled by 
a factor I (i.e. IA = {x : ^ _1 x 6 A}). In this limit semi- 
classical methods are often useful. To study the eigenval- 
ues of Ke,ia, it is convenient to characterize the behavior 
of Tr/(ifg ia) for various functions /. We immediately 
see, since in our situation [Hq(p),A] is supported only 
on the boundary of the system, that in the volume term 
we may commute the \A terms in the perturbations W* 
and the free propagator. This gives the asymptotic lead- 
ing term for the number of sub-gap states as the number 
of eigenvalues of Kejo, larger than one. The semiclas- 
sical approximation expresses it as an integral over the 
classical phase space: 

N(E, 0) - l d vol(A) [&[(H (p) - E)- l a x ) - l]d d p(24) 



where is the Heaviside step function. 

We may view the operator Ke,a as a higher dimen- 
sional block of a Toeplitz operator: the kernel depends 
only on x — x 1 , however the x, x' indices are restricted to 
be in A. This is equivalent to chopping a "block" from 
a regular matrix. The theory describing the asymptotic 
properties of blocks of Toeplitz matrices is based on var- 
ious Szego theorems, and appears in numerous problems 



in physics. Most famously, the strong (two-term) Szego 
limit theorem was initially used in the celebrated compu- 
tation of the spontaneous magnetization for the 2D Ising 
model by Onsager (see e.g. 23J). The Szego limit the- 
orem also plays a special role in entanglement entropy 
studies, see in particular [24| [25]. 

To try and describe the next term, we turn to a higher 
dimensional generalization of the Szego limit theorem, 
which has been extensively studied by Widom [55]. De- 
note Ge(p) = (Hq(p) — E)~ x o x . Given a unit vector di- 
rection n x , we can define an operator acting on functions 
of a single variable, with kernel : 



Qe,: 



t) 



,ip|l(s-t) 



Ge(p± +n ar p||)dp||. (25) 



Using [55] we see that 

Trf(K EM ) ~ ^§#1 / f[g E (p)]d d P + 
(2^r J aA dxf p± , nx=Q d d - 1 p ± tr[f{ X+ gE,n !C ,p ± X+) 

-X+f(GE,n^ P Jx+}+o(l d - 1 ) (26) 

where n x is the normal to dA at x. For the case of 
an s-wave superconductor the operator QE,n x ,p ± can be 
computed explicitly: 

Ge,ti x , p A s — *) = 



2(Ag-V) Im ( k + ( P± ) )( A o + **■*) 

«|fc + (j>X) N 



+i-J%A Re('- 



*+(pj-) 



E 2 , and 



where where k±(p±) = yp|-fi±i\/Ao' 
Refc ± >0 (here ia y = (0,1,-1,0)). 

In our case, we are interested in the number of eigen- 
values above 1. Thus we need to solve this problem for 
/(A) = 9(A — 1). However, for this function the argu- 
ments leading to (26 1 arc no longer applicable, but rather, 



due to the discontinuous nature of /, the area term as 



written in (26) is divergent. We conjecture that in such 
a case, the boundary term should be modified to: 



id-l. 



W)^ fdA dx / Pi .n«=0 d ° >- ,] " 



-xi+f{Gi 



.)xm 



f(Xl + GE,n x ,p ± Xl+) 

-o^- 1 ), (28) 



where xi+( s ) = 9(0 < s < I). In analogy to the 
appearance of logarithmic corrections when consider- 
ing similar expressions for the entanglement entropy of 
fermions in d-dimensions, we expect this term to scale 
as 0(l d ~ 1 log?), which also seems to be consistent with a 
preliminary numerical investigation. 

In this case, (as well as for any isotropic superconduc- 
tor) , this integral is independent of the direction n x , and 
so we have that the boundary term is given by: 



L^l V ol(dA)m f (E) 



(29) 



for some function rrif(E). The transformation / — ► wif is 
independent of the shape, and we believe can be studied 
using Wienner Hopf methods for specific cases. 

Finally, we remark that it has been argued that there 
is a mini-gap present in the local density of states if 
the shape of the normal region corresponds to a chaotic 
"Andrcev billiard" [27l430| . The mini-gap is of order 
Ao \og(kpl), and it's nature depends on te/te, where td 
is dwell time in the normal region and te is the Ehren- 
fest time. This question has recently been revisited in 
|31j . To arrive at this result, the strategy is to consider 
the BdG equation within the Andreev approximation and 
analyzing possible classical trajectories. It would be of 
interest to see whether such a result may come out of a 
Birman-Schwinger like analysis. 



VI. CONCLUSIONS AND OUTLOOK 

To summarize, in this paper, we studied the bound 
states in superconductors using a Birman-Schwinger ap- 
proach. The advantage of the method is its validity for 
general perturbations, and beyond Andreev approxima- 
tion. Our main results are stated in following equations: 
The bound on sub-gap states (11), the Cwikel-Lieb- 



Rozenblum type inequality for BdG ( 20 1 and the asymp- 



totic expression for the number of states in large regions 
with the conjecture (28) for the scaling I 1 \og(kpl) of 
the sub-leading term. For most physical situations the 



Birman-Schiwngcr bound ( 11 ) seems a better bound than 



(2'J the CLR bound (20), since it gives a lower upper bound 



for large fermi energy \i. 

In the same way it is possible to study the effect of 
adding disorder potentials, as well as spatially varying or 
momentum dependent A. Because of the wide interest 
in the properties of vortices in such systems it is also 
of great interest to understand if the Birman-Schwinger 
methods can be useful in the study of bound states in 
vortices. Indeed, it is not immediately clear how this 
should be done: The presence of core states in vortices 
seems to be mainly due to the phase winding around the 
vortex core (see, e.g. [32]), and arguably less sensitive 
to the supression of order parameter in the core. Since 
the phase winding is of a topological nature (i.e. it 
is discrete) it cannot be considered as a perturbation 
which is turned on in a continuous way from zero, as 
in the case we have considered. The present approach 
may be utilized to understand the interplay between the 
topological and non-topological state binding. To do so, 
one should take as the unperturbed vortex an idealized 
radially symmetric vortex, and then add phase gradients 
and order parameter suppression as perturbations, and 
study the additional bound states introduced in this 
manner. 
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d/2 



Let Vd — r(d/2+i) 
sphere. Then: 



be the volume of a d dimensional 



Appendix: Cwikel-Lieb-Rozenblum for BdG 

A. Estimates on singular values 

To control the spectrum of the Birman Schwinger oper- 
ators, we will use the following inequality (p > 2), (Thm 
XI.22 in [33]): 



9n(tA)<n-^ p C Ptd \\a\\ 



p\\v\\p,w> 



(30) 



where s n is the n-th singular value of the operator 
ab = a(x)b(iW), i.e. the nth largest eignvalue of the 



operator y^ab^ab. Here C p ,d is a constant, ||a|| p = 

(Jdx|a(x)| p ) 1 / p is the LP norm, and the value ||6|| PiU , 
is defined as follows (pg. 30 in [34]). Defining: 



M b (t) = vol(x\\b(x)\ >t) 

\ P ,w is given by: 

C p 



\ PtW = inf{C : M b (t) < —Vt} 



(31) 



(32) 



Note 



,. w is not strictly speaking a norm, since it does 
not satisfy a triangle inequality. It is called weak since 
L p C L™. 

To familiarize ourselves with computing ||6|| PiU , con- 
sider first the following example: 

Example: Let us compute ||&||3,«, for b(x) — (x 2 + 
c) -1 / 2 assuming c > 0. We have: 

M b (t) = V0l(x\(x 2 + C)- 1 ^ > t ) = 

vol{x 2 < £ - c)6{^ - c> 0) (33) 

In 3d this gives: 



M b (t) = f\^-c\^ 2 i 



:?-c>°)<y* ^ 



since the above inequality approaches equality for t — > 0, 

3 



the infimum in (32), is given by 



3,iu 



(f) 1 ' 3 . Note 
that ||&||3,to does not depend on c. 

In the derivation of the Cwikel-Lieb-Rozenblum in- 
equality we will need the | |fo| |3,iu estimate for 

b(iV) = ((-V 2 -m) 2 + cy 1/2 , (35) 

with m, c > 0. Here we show that 

\\b\\ 3 . w < c-Hv 3 )^max\l, ( 324 !^) 1 / 12 j (36) 

L v c ' J 

As in the example above, we define 

M b (t) = vol(x\{{x 2 - to) 2 + c)- 1 ' 2 > t) 

= vol((x 2 - m) 2 < i - c)9(± - c> 0) (37) 



M b {t) 



'v d (m+^-c) d l 2 



if 



Vd 



(to + 



c) d ' 2 - (r 



< t < (to 2 + c) _1/2 

~) d / 2 ] if 
(to 2 + c)" 1 / 2 << < 4 



v^ 



if 



t > 



(38) 



From the small £ behavior, we see that to bound M b 
uniformly as in the definition ( |32| the power has to be 
at least l/t d / 2 , i.e. (p > d/2). In the first sector we see 
also that the constant has to be at least Vd- However i>d 
may be not good enough since the bound in the region 2 
may be different. Let us estimate it: 
For (to 2 + c)- 1 / 2 <£< X, 



t'd 



z y/2 _ ^ 



(to + 

v d j\dsd s ( m + s^-c) d / 2 = 



)d/2 



2 



y / £-cf\ds(m + S y , £-c) d / 2 - 1 < d ^(2m) d / 2 - 1 



where we used that in this region to > « / -h — c. 

Now we have to write this as bounded by a l/t d ' 2 
expression. If d/2 < 1 we are done, but if not, d/2— 1 > 0, 
like d = 3 case, using -K= > 1 we have: 



i 



±{2m) d / 2 ^ < d ^(2m) d / 2 - \ d/2 _ 1{ y c)d/2 _ 1 = 

dv d (2m,) d ^- 1 
(v^) d / 2 - 1 t d / 2 

Combining all these results we see that (for d > 2) 



(39) 



M b (t) < max(l, 



d(2m) d / 2 -\ v d 



{yfcY' 2 



r) 



t 3/2 



-3/2 



And thus we see that in this case we are in the weak L v 

space, i.e.: ||6|| 3 /2,iu < (v d max(l, 3(2 ™) 4 )) 2/3 . 

The Cwikel estimates require p > 2. Let us estimate 
||6||3 ]U) . Since t~ 3 / 2 < i -3 , say, for t < 1, we need to deal 
with the large t behavior of M b . Let us use the following: 



Lemma: If M b < Ct 



and M b = for t > to, then we 

(40) 



have, for any ct\ > a that: 

||&IU,™<(C~ Q c) 1/p - 

Proof: for t < t we have 

M b (t) < cr a < ct- a (t a /t) ai - a < ct- ai Q- a 



and M b = for t > t therefore M b < Ct%~ a H~ ai for all 

t > o. n 



In our case Mb drops to zero when t = c 1 / 2 . There- 
fore, for b(x) = {{x 2 — to) 2 + c) -1 / 2 we have, for p > 3/2: 

H&IU < c-^cl(v d ma X (l, 3{2 ^ /2 )y/i> (41) 
and in particular: 



i/ si/t l\ ,324m% 
|3,td < c ^(v 3 ) 1/s max 1, ( ) 



, 324m 2 N i/i 
c 



(42) 



B. Derivation of the main inequality 



Here we repeat almost exactly the analysis of Cance- 
ller, Levy-Bruhl and Nourrigat, "remarks on the spec- 
trum of dirac operators" [25]. The main difference be- 
tween our work and theirs is the presence of the chemical 



potential, which appears in the estimates such as ( 60 ) 
Assume ||V - V A || < |A |A/4, then 



\\(H + V)u\\ < A VT^X\\u\\ ^ 
||(iJo + V*H<A OA /l-A/2||u|| 

To check this note that: 

\\(Ho + V x )u\\ 2 < 

\\(h + v)u\\ 2 + 2\\(h + v) u \\\\(Vx-v)u\\ 
\\(V-V x )u\\ 2 < 



(43) 



|A | 2 (l-A) + |A | 2 |yr^A+|A 



12 Al" 

I 16' 



<|A | 2 (l-A+f % /T^A+^)|M| 2 < 



|A | 2 [1 



HP 



(44) 



The last inequality can be checked by differentiation. 
Thus we have: 

||(F + ^) U || 2 -|A | 2 [1-|]|H| 2 <0^ 
(u\(H 2 - |A | 2 ) + [H , V X } + + V 2 + |A | 2 f )|u) < 

Noting that for BdG H 2 = (-V 2 -/x) 2 l + |A | 2 l, we get 
for such us: 



H(-v 2 - M ) 2 



|A | 2 H5a|^<0 



(45) 



where S\ = [H n , V\]+ + V 2 . Therefore we are interested 
in the number of bound states of (—V 2 - \i) 2 + |A | 2 f + 
S\. These are given by a Briman-Schwinger method. 
Observing the Birman Schwinger kernel: 



K 



BdG 



[(_ V 2-^)2+|A pA]l/= S 'A [( _ V 2_ /j) 2 + | Ao |2A]l/ 2 (46) 

We now note that the S\ contains several terms, of the 
following forms (there are also matrix indices, but let 
us ignore these, since we are going to add the estimate 
for each element). We also assume A real for simplic- 
ity, also, here the perturbation is actually of the form 



-Ai(a;)<7 x , so we appropriately define A A as we defined 



V x 



Next we use the Cwikel inequalities ( 30 ) for p — 3 and 



the standard inequalities for singular numbers: 
Fan's Inequality: 



and 



[X n+m +i(A + B) < Hn+l(A) + n m+ i{B) 



Hn(AB) < \\B\\n n (A) 

fin+m+l(AB) < fl n+1 (A)fl m+1 (B) 



(47) 



(48) 



to deal with the different terms in KbcIG- 

We get estimates for the n th singular number of three 
different forms (here b is as above in (35)): 

1) 



Ai = A a 



[(_ V 2-p)2+|Ao| 2 |] 1 /2 V Aa(«) X 



V AaW [ (_ V 2_ ai) 2 + | Ao |2 AJ1/2 



*n(Ai)<Ao(||VAA|| 8 C3,d||6||8,. 



(n/2)!/s 



(49) 



(50) 



Here we used standard inequalities for the singular val- 
ues of sums and products of operators like (48), to have 
s n {AB) < s [n/2] (A)s [n/2 ](B). 
2) 



A 2 - A 0[( _ v2 _ M) 4 |Ao|2 x ] i /2 A A (a:) x 

Gives in the same way: 

1 



s n (A 2 ) < (|| 



AxhCsAblkw 



(n/2)Vs 



(51) 



(52) 



3) 



A, = i 

> 3 [(-V 2 - M ) 2 + |A | 2 A]i/ 



^(-V 2 -/i)x 

Aa(z) [(_v 2 -m) 2 + |A |2|]i/2 



(53) 



l IF ^ 2 -^___(_v 2 - //)|| < 1, we have: 



s»(40<||A A || 3 c M ||&|k. 



(n)V3 



(54) 



Combining these estimates and adding them, (and noting 
that ||v / Ax|| 2 = ||A A || 3 / 2 ) we conclude the inequality: 

Sn {K BdG ) < C"[Ao||A A || 3/2 ||&|| 2 !U; -i I + 

l|A A ||l||6|| 2 ^43+H A >ll3ll fo ll3,-^l ( 55 ) 
' n i s n ' . 



Now we can get the inequality as follows: we are looking 
for the largest N for which Sff > 1. Thus we multiply by 
TV 2 / 3 the above inequality: 



N 2/3 < C 



AoIIAaIM 



\3,w 



IAaIIII 



IA 



A3 



11,™ 4 

,7V 1 / 3 



(56) 

(57) 



Thus if N m satisfies: 



N% 3 = C 


A ||A A || 3/2 ||6||^ + 


±x\\l\\b\\l w 


+ ||A A || 3 ||6||3, ro iV 1 / 3 



Number of bound states must have N < N m . Let us 

2/3 1/3 

bound N m from above. Writing 7V m = A + BNm =>■ 
N i/3 = B+VW+a _ Since ^ B 2 + AA < b + 2\/J, we 

have: Nm < B + y/A. Furthermore, using the general- 
ized mean inequality (|o| + |6| + \c\) k < 3* _1 (|a|* + |6| fe + 
\c\ k ) we have: 



N m < C 



(Ao||A A || 3/2 ||6||^) 3 / 2 + 



|Aa|||| 



C"\ 



li,J 3/2 



IA 



A3 



\Z,w 



A 3 / 2 || A, || 3 / 2 



h.wY 



*A||s 



< 



(59) 



In the BdG context we have m = \x and c = A|A | 2 in 



( 35 1 . There are two cases: 



1) fi > v ° p- — : We have we have 



|s,« < 2( V3 ) 1/3 ( 



/3( 9 M xl/6 



2A 4 |A | 



Using ( 60 1 we have 



IL<M; 9 " ^ 



^2A 4 |A | 



(58) So finally: 



,,1/2 

M < C—tL—— 

JV m V O A 2| A() |2 



v 3/2 



1 3/2 



AS"||A A ||^ + ||A A ||g (62) 



for some C 

We can now write this explicitly as an integral: 



N m < C 



/' 



1/2 



A 2 |An| 2 



d.i' 



OP J\5A{x)\>^>t 



.3/2 



3/2 



(60) 



(61) 



A° "(6A x y^ + (6A X ) 



The optimal constants in ( 30 1 are not known in general, 



but the various upper-bounds for them are all of order 
one. Finally, relating A and E we get (20 1. 



/ a 2_ o2 

2) Similarly, for [i < v ",- — and using (42 1 we get 



(211 
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